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The mean-first-passage-times (MFPTs) for a vacancy that diffuses 

(via one- and six-jump cycles) in a two dimensional ordered binary 

alloy are evaluated using the properties of random walks on networks. 

cn . We investigate the effect of temperature and relative barrier height 

on the ratio between the MFPTs of the two cycles. At low tempera- 
ture we find that the six-jump cycle takes shorter time while at high 



^^ ■ temperature the one-jump cycle takes shorter time than that of the 

OO , six-jump cycle for the range of parameters considered. 

O 

1 Introduction 

H, 

The mechanism by which a single vacancy diffuses in mono-atomic 

crystalline material is basically through site exchange with one of its 

nearest-neighbor atoms. The vacancy diffusion continues through the 

material successively in a random way such that each new site occupied 

is usually energetically identical to any other earlier occupied site. As 

such, order is maintained throughout the vacancy diffusion. 
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Figure 1: A binary alloy on two-dimensional lattice with a single vacancy 
occupying a site on sublattice A. 



The situation is not so simple if the crystalline material is, for 
instance, composed of a binary alloy which consists of two interpen- 
etrating simple cubic sublattices that are predominantly occupied by 
two different atoms, A and B. Figure 1 is an illustration of a binary 
alloy in two-dimensional lattice with a single vacancy. Whenever the 
vacancy on site A exchanges its site with one of its nearest-neighbor 
atoms, the process leads to disorder in the crystalline structure. If the 
vacancy randomly moves successively via nearest-neighbor jumps, a 
string of anti-structure atoms would lead to disorder in the material. 
To avoid this problem of disordering, two alternative diffusion mecha- 
nisms are likely to take place: either jumps of the vacancy to further 
distant sites on the same sublattice or a cycle of successive intermedi- 
ate jumps to nearest-neighbor in which the atomic disorder appearing 
during the earlier part of one cycle is followed by successive healing 
during the later part of the cycle. The first alternative is usually called 
one-jump cycle. The prominent cycle for the second alternative is the 
six-jump cycle. It was first suggested by H. B. Huntington and later 
discussed by Elcock and McCombie pQ. Figures 2a and 2b illustrate 
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Figure 2: A possible path for (a) six-jump cycle and (b) one-jump cycle. 



the possible paths for six- and one-jump cycles, respectively, for a 
binary alloy in a two-dimensional lattice. 

A previous work dealt with calculating diffusion coefficient via six- 
jump cycles using rate equation method [2]. An even earlier work 
[3] used mean-first passage method to compare the diffusion coefficients 
of the A and B atoms. In our present work, we raise a specific ques- 
tion regarding the times required for a vacancy to diffuse via the two 
dominant cycles and compare them using appropriate parameters. 

Which one of the two alternatives will the vacancy prefer to diffuse 
through the lattice? One way to answer this question is to compare 
the times taken for the vacancy to go from one stable state to the next 
stable state via the two alternatives. 

The problem of vacancy jump from one site to the next through 
site exchange with an atom can be seen as a barrier crossing problem 
of an idealized Brownian particle representing the combined vacancy- 
atom pair involved in the site exchange. The direction of vacancy 
motion can be taken as the direction of movement of the Brownian 
particle. Since each state to which the vacancy jumps is either stable 
or metastable, it is reasonable to consider the next jump process as 
independent of its previous one. In other words, we will consider the 
Brownian particle to have enough time to get thermalized at each site 
before taking the next jump. 

The problem of evaluating the times the vacancy takes from one 




Figure 3: Schematic diagram showing energy profile for (a) six-jump cycle 
and (b) one-jump cycle. Note local jump probabilities p (up) and q (down) 
in (a) and jump probability po in (b)). 



stable site to the next stable site via the one- and six-jump cycles de- 
pend on the corresponding potential energy profiles of the two alter- 
native paths. These times are usually called mean- first-passage-times 
(MFPTs). In the case of the one-jump cycle, the vacancy has to make 
a single jump over a high barrier to the nearest site on the same sub- 
lattice. In the case of the six-jump cycle, the vacancy has to undergo 
through six local jumps where the total barrier height is subdivided 
to ultimately arrive at a nearest site on the same sublattice through a 
longer path than that of the one-jump cycle. Figs. 3a and 3b are the 
energy profiles for the six and one-jump cycles, respectively, during 
the course of each cycle. Given the specific potential profile for the 
one-jump cycle, one can use the standard method of solving Brownian 
diffusion in a potential field to evaluate the MFPT for the one-jump 
cycle. On the other hand, to evaluate the MFPT for the six-jump 
cycle we use a technique that has first been formulated by Goldhirsch 



and Gefen [H [5] and later applied by one of us in collaboration with 
others [6j. This technique requires knowledge of local jump proba- 
bilities between successive sites as an input in order to evaluate its 
MFPT. 

In this work we consider the vacancy diffusion on a two-dimensional 
lattice. In this case, there are four possible nearest stable sites on the 
same sublattice to which the vacancy can jump by either performing 
the one- or the six-jump cycle. In principle, the vacancy can attempt 
to diffuse via all these paths, select one of the paths and ultimately 
reach one of the four sites in one cycle. 

The rest of the paper is organized as follows. In section 2, we con- 
sider the multiple paths scenario for the two alternative jump cycles 
that describe vacancy diffusion on two-dimensional lattice and deter- 
mine their corresponding closed form expressions for their MFPTs. 
Section 3 compares the values of MFPTs as a function of some of the 
parameters of interest and explores the different possibilities defining 
them. In section 4, we present the summary and conclusion. 

2 MFPTs along the one- and six-jump 
cycles 

The MFPT,n, that the Brownian particle takes to complete the one- 
jump cycle along a single path is worked out in the Appendix. We 
consider a potential profile similar to Fig. 3b but for simplicity the 
potential profile is considered to be a piece-wise linear potential as in 
Fig. 6 of the Appendix. The expression for the MFPT in the low 
temperature regime turns out to be (see Eq. (A10) in the Appendix) 

T1 = D \~W) eB (1) 

where a is the distance between nearest neighbors, Eq is the barrier 
height energy, D is the diffusion coefficient, k B is Boltzman's constant 
and T is the temperature of the crystal medium. One should note 
that there are four possible stable lattice sites to which the vacancy 
jumps on the same sublattice. If the vacancy chooses the one-jump 
cycle, the MFPT T14 it takes to reach one of the four possible sites on 
the same sublattice is simply 

n 1 (k B Ta\ 2 *%, 

ri4 = T = 4iH^rJ e ■ (2) 



On the other hand, if the vacancy chooses the six-jump cycle to ul- 
timately arrive at one of the four stable lattice sites on the same 
sublattice there are two possible paths. Therefore, there are a total 
of eight possible paths which the vacancy can select to move by one 
sublattice distance. The evaluation of the MFPT, tqq, the vacancy 
takes to move by one sublattice distance via the six-jump cycle after 
choosing one of the eight possible paths is done using the technique 
formulated by Goldhirch and Geffen [U [5] and given by 

-q{p + q) f2p 3 -7p 2 -9pq + q 2 \ 
T68 4p V<7 5 q 2 {p + q) 2 )' l ' 

where p and q are the local jump probabilities (up and down) over the 
rugged potential corrseponding to the six-jump cycle (see p and q in 
Fig. 3a) The closed form expressions for p and q are derived in the 
Appendix. 



3 Result and Discussion 

The expressions for the MFPTs are functions of a few physical quan- 
tities and include the barrier heights of the one- and six- jump cycle, 
Eq and E\, the background thermal energy, ksT, of the binary al- 
loy, the local barrier heights, e± and 62, for the six-jump cycle (see 
Fig. 6 in Appendix), the diffusion coefficient D and the spacing, a, 
between the sublattices. In order to compare the MFPTs between the 
two alternatives we first identify two dimensionless quantities that are 
parameters controlling the MFPTs. One is the dimensionless quantity 
b = y±j; that controls the local hopping rate for the six-jump cycle. 

The other dimensionless quantity is the ratio r = -^- which compares 
the total barrier corresponding to the one-jump cycle with that of the 
total barrier height for the six-jump cycle (see Figs. 3 (a) and (b)). 

Let us define a dimensionless quantity, f m , that compares the 
MFPT, T14, for the one-jump cycle to that of the MFPT, tqs for the 
six-jump cycle given by 

/ m = |log(— )|- (4) 

T68 

Fig. 4 shows three plots of f m as a function of b (scaled local bar- 
rier height for the six-jump cycle) corresponding to three different r 
values that compare the total barrier heights of the two cycles :r = 
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Figure 4: Plot which shows f m versus b for fixed r = 0.8, r = 1 and r = 1.2. 



0.8, laradl.2. Each plot has a region of negative slope as well as a 
region of positive slope. Note that the region of negative slope corre- 
sponds to a situation where the MFPT, r 14 , via the one-jump cycle 
is smaller than that of the MFPT, tqs, via the six-jump cycle. This 
implies that in the high temperature regime the one-jump cycle takes 
shorter time than that taken by the six-jump cycle. This is because 
when the background thermal kick is high enough, the vacancy can 
easily cross the barrier height through one-jump cycle. On the other 
hand, the region of positive slope corresponds to a situation where 
tqs is smaller than T14. This implies that at low temperature the 
six-jump cycle takes shorter time than that of the one-jump cycle. 
This is because, at low temperature, the background thermal kick is 
weak for the vacancy to cross the barrier via one-jump cycle. On the 
other hand, for the six-jump cycle there are six small barrier heights 
which can be crossed with relatively small thermal kicks. Thus at low 
temperature the vacancy can cross the relatively small local barriers 
quickly compared to the large barrier of the one-jump cycle. At the 
inflection point of the plots, f m is zero and this corresponds to a sit- 
uation where the MFPT via both cycles is the same. This clearly 
shows that at a certain temperature the MFPT for a vacancy diffus- 
ing through one-jump cycle is equal to that of vacancy diffusion via 
the six-jump cycle. Comparing the three plots, when r = Eq/E\ gets 
large the MFPT via the six-jump cycle predominantly gets relatively 
shorter compared to that of the one-jump cycle. Fixing 6 = 2 and 
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Figure 5: Plot which demonstrates f m versus r for fixed b = 2.5 and b = 2. 

b = 2.5, we investigate how the f m behaves as a function of r as 
shown in Fig. 5. The figure clearly demonstrates that the MFPT 
for one-jump cycle becomes shorter than the MFPT of six-jump cycle 
when the thermal background temperature is strong enough as long 
as the potential barrier for one-jump cycle, Eq, is smaller than E\ for 
six-jump cycle. 



4 Summary and conclusion 

We studied the MFPT for a vacancy diffusing in two dimensional 
binary alloys. We investigated how the MFPT behaves as a function of 
the two model parameters. The central result we obtained shows that 
the six-jump cycle takes invariably shorter time than that of the one- 
jump cycle for the parameter ranges considered at low temperature 
regimes. When the back ground temperature is strong enough, the 
one-jump cycle is dominant. 

As a concluding remark, there are two issues we would like to point 
out for future consideration. The first one concerns the need to explore 
other ranges of parameters not considered in this work. The second 
one is the need to relate our findings to experimental results. 
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Appendix 



In this Appendix we derive the values for p, q and po in terms of the 
model parameters. For a Brownian particle that moves in a highly 
viscous medium under the influence of external potential V(x), the 
Langevin equation that governs the dynamics of such a particle is 
given by: 

where 7 and T denote the constant friction coefficient and the tem- 
perature, respectively. V(x) is the external potential, ks is the Boltz- 
man's constant and £(£) is the delta correlated noise term. The cor- 
responding Fokker-Planck equation can be written as, 

d t p(x, t) = d x [^lp(x, t)] + Dd 2 x p(x, t). (A2) 

7 

Here p(x, t) represents the probability of a particle to be found at 
a position x at time t, and D = ksT/j is the diffusion coefficient. 
Starting from Eq. A2, one can derive the expression for MFPT in a 
bistable potential whose inverse is the jump probability [7]. Taking a 
piecewise linear potential which is described by 



V(x) 



2a ' x — U ' 

3gig < X < — • 

4|k + 2ei , f< x ~<% (A3) 

^f-ei, x>a, 



and noting that €2 = ei/2, the MFPT r(0 — ► a) to reach x 
starting from x = is 



-n —)(A—Ae h B T +8e 2k B T — e 2k B T ( 

9Dej A y k B T 

(A4) 



r(0 — > a) = ( nrk2 )(4-4e fc BT+8e 2fc B^-e 2fc B^( r ^ ; +8)). 
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Figure 6: Asymmetric piecewise linear double well potential. 
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When ei is large compared to k B T, the above expression takes a simple 
form: 

4 / k B Ta \ 2 _a_ , A N 

r(0 — a ) = - M— e k * T . A5 

For the case ei is large compared to k B T, probability p to jump from 
x = to x = a is the inverse of r(0 — > a) so that 

'-t(i5t)''*- < A6 » 

Similarly for the reverse case, the MFPT r(a — ► 0) for the vacancy 
to reach x = from x = a is given by 

2a 2 e k B T (kuT) 2 _!2_ _!£2_ 2ia_ 2e 9 

r( Q ^0) = ( Za 9jD ^ 2 J )(4-4e^+8e^ +e ^(.i!|- 8 )). 

(A7) 
For low temprature regime the above equation takes a simple form: 

4 / k B Ta \ 2 _a_ , 4 N 

D V e i + e 2/ 
and the associated local jump probability q is given by 

If the bistable potential is symmetric with large barrier height Eq 
(compared to thermal energy) and width a, the MFPT taken to jump 
in both directions is the same and is given by 

Tl = ( (0^a) = i(^) 2 e * (A10) 

while the corresponding jump the probability, Po = q = p takes the 
value 
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